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0\ '. Abstract 

o> ■ . . 

A simphfied Lagrangean closure for the Navier-Stokes equation is used to study the production of inter- 
s' ^ mittency in the inertial range of three dimensional turbulence. This is done using localized wavepackets 
Q ■ following the fluid rather than a standard Fourier basis. In this formulation, the equation for the energy 
transfer acquires a noise term coming from the fluctuations in the energy content of the different wavepack- 
ets. Assuming smallness of the intermittency correction to scaling allows the adoption of a quasi-gaussian 
^*| I approximation for the velocity field, provided a cutoff on small scales is imposed and a finite region of space 
is considered. In this approximations, the amplitude of the local energy transfer fluctuations, can be calcu- 
lated self consistently in the model. Definite predictions on anomalous scaling are obtained in terms of the 
^ ■ modified structure functions: << E{l,a) >%_>, where < E{l,a,r,t) >j, is the part of the turbulent energy 
OO I coming from Fourier components in a band (a — l)fc around k ^ l~^, spatially averaged over a volume of 
size R ^ — ^- around r. 
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I. Introduction 

The statistics of large Reynolds numbers 3D (three dimensional) turbulence is characterized by scaling 
behaviors of the structure functions: S{1, q) =< vf >, Vi(r, t) = v(r + 1, t) — v(r, t), which, in first approxi- 
mation, appear to follow the Kolmogorov relation: S{l,q) =(x P', Q ~ q/3 [1]. However, both experiments 
[2] and numerical simulations [3] show the presence of corrections to Kolmogorov scaling, which become 
more and more pronounced for higher order moments. These corrections go in a direction of increasing 
nongaussianity of vi as I ^ 0: S{1, q) ex l^" , with S(q = Cq — Cqi satisfying the convexity condition: ../ < 0. 
Of course the relations above hold only for scales corresponding to the inertial range (and a part of the 
viscous range [4] if only ratios of moments are considered), so that for finite Reynolds numbers: Re < oo, the 
generalized kurtosis K{2q) = [S{l,q)/ S{l,2y]i^„ is finite. However, a situation with scaling corrections in 
the form just described, persisting for Re -^ oo, would imply that K(q > 2) ^ oo in this limit, corresponding 
to infinitely intermittent small scale velocity fluctuations. 

The traditional approach to the issue of intermittency dates back to the remarks by Landau [5] on 
the presence of spatial fluctuations in the energy transfer between scales and its effect on the turbulent 
dynamics. From the Refined Similarity hypothesis of Kolmogorov [6], down the line to the Beta Model 
[7-8], the main ingredient is the assumption of a turbulent dynamics acting locally in scale, so that the 
energy density at a given point in space is essentially the product of coefficients, each with an independent 
fluctuating component, and each describing the transfer of energy between successive, contiguous scales. In 
this picture the source of intermittency lies in the inertial range and acts locally in scale; therefore, the 
scaling corrections is predicted to be independent of Re. This mechanism can be interpreted as a rule for 
the construction of a multifractal [9-10] in which the measure is given by: (i/i;(r, i) ~ I < (dxVx)'^ >i (r,i), 
where < >/ (r, t) indicates a spatial average taken at time t in a an interval of length I along the x-axis, 
centered around r. In this way, the multifractal dimensions Dq, defined by the relation < d/x^ >~ l(<}-i)Dg^ 
are given by: {q - 1)(1 - Dq) = qC2 - C2q- 

Recently, an alternative explanation for the presence of scaling corrections has been proposed, namely 
that what is observed is a finite size, i.e. a finite Reynolds number, effect [11-12]. This point of view is 
supported both by the smallness of the corrections and the strongly nonlocal character of finite size effects 
with respect to scale, as shown in numerical simulations of both Navier-Stokes dynamics [2] and reduced 
wavevector models [11]. In [12] a dynamical explanation for this effect was presented, based on closure 
arguments. The basic idea is that the strong intermittency of the viscous range, due to the interplay 
between dissipation and nonlinearity in the dynamics of small vortices [13], may be enough to generate, in 
the inertial range, scaling corrections of the appropriate size for all Reynolds numbers of practical interest. 
In this picture, the total intermittency, as parametrized by K{q), would remain finite, while the scaling 
corrections would tend to zero in the infinite Reynolds number limit. 

In this situation, it would be interesting to have some quantitative assessment of the magnitude of 
the finite size effects, relative to the amount of intermittency produced in the inertial range by energy 
transfer fluctuations. However, while [11] and [12] can afford some quantitative prediction on the size of the 
first effect, all models dealing with the second are heavily phenomenological and have no connection with 
real Navier-Stokes dynamics. An exception is the very recent paper by Yakhot [14], in which quantitative 
predictions have been obtained after approximating the turbulent energy dynamics with that of a passive 
scalar. 

The purpose of the present paper is to obtain quantitative predictions on the i?e-indepcndent part of 
anomalous scaling (if present at all), following the route of evaluating the energy transfer fluctuations from 
statistical closure of the Navier-Stokes equation. The basic difficulty in using this technique is the real space 
nature of the quantities that have to be calculated: the moments S{1, q) [15], while the main component of 
the turbulent dynamics, the energy transfer towards smaller scales, is better described in a Fourier basis. 
Although a real space closure of the Navier-Stokes equation was derived in [16], we found very difficult to 
extend the model to study energy fiuctuations, especially when trying to separate contributions at different 
scales. 

Here we prefer to follow the idea of Nakano [17] and, in a different context, of Eggers and Grossmann 
[18], of using a localized wavepacket representation of the Fourier- Weierstrass type rather than a global 
Fourier basis or a purely real space one. In this way, energy fiuctuations on a scale i?, instead of being 



buried in the phase relationship between Fourier modes separated by Afc ^ R^^, are described by the space 
dependence of the wavepackets. Now, if the wavepackets are centered around wavevectors lying in shells of 
radii A:„ = a^ko, the spatial extension of a wavepacket at fc ~ l~^ will be R{1) ~ Z/(a — 1). This means that 
a new object is being used in place of S{l,q) to give a measure of intermittency: the generalized structure 
function S(l,a,2q) —<< E{l,a) >%>, where E{l,a,) = E{l,a;r,t) is the total energy of Fourier modes in 
the shell /c„ ^ l^^, and R — R(l, a) — l/{a— 1). More in general we shall consider the situation in which also 
the quantity c^ = / _-^, „ is treated as a free parameter, resulting in the definition of still another structure 
function S{l,a,Cw,q) interpolating between the intermittency free limit S{l,a,0,q) and the original case: 
Sil,a,cT''.q) = S{l,a,q). 

The shell width (a — l)fc plays here a crucial role; this can be understood better by looking at how 
energy is transferred between shells as a rule for the construction of a multifractal. In this picture, a small 
value for a — 1 corresponds to multipliers between one scale I and the next l/a being constant on domains 
R{1) > I (in typical examples like the various "middle third" Cantor sets [10], one has R{1) ~ /, i.e. the 
fluctuation scale and that of the geometrical structures is the same). In terms of turbulent dynamics, this 
has the following interpretation: R(l) is the scale of the eddies which contribute the most to the straining of 
those of size I; the parameter l^^R{l) = (a — 1)^^ gives therefore the degree of nonlocality of the nolinear 
interaction. In principle there will be fluctuations in the transfer of energy between eddies of size / also 
at scale smaller than r{l) and this effect will result in correlations between phases of different wavepackets 
separated by Afc > R^^. This effect contributes to the scaling of S{l,q), but not to that of S{l,a,q), and 
the last one is likely to underestimate the actual value of SCq ■ 

Now, recent analysis performed by Domaradzki and Rogallo [19] on numerical simulations and statistical 
closures of the Quasi Normal Markovian type [20] has shown that there is indeed a separation of scales 
between straining flow and strained eddies, so that fluctuations in the transfer, occurring at scales similar to 
those of the modes exchanging energy, are not expected to be large, implying: S{1, q) ~ S{1, a, q) even for a 
close to 1. More importantly, it is this separation of scales that allows to consider meaningful wavepackets, 
extending over several characteristic wavelengths, rather than having to deal with the usual eddy breaking 
picture, that is very attractive on grounds of simplicity, but does not allow to make any contact with the 
Navier-Stokes dynamics. 

There is a second conceptual difficulty in using closures to study multifractal intermittency. This is 
the contradiction between an ansatz of quasi gaussianity and the "infinitely non gaussian" character of a 
multifractal with no cutoffs, as attested by the equation: K{q > 2) = oo. A quasi gaussian hypothesis 
becomes meaningful however when studying limited regions of space and limited ranges of scales, which 
is possible if the degree of nonlocality of the nonlinear interactions is not too high and the intermittency 
correction S(^q is sufficiently small. In this way, although the distribution of values of, say dxVx, over 
different averaging volumes Vi is infinitely intermittent (if the total volume Vtot ~ Ul^ is infinite itself), 
the moments of dxVx in each Vi will be close to gaussian [if the ultraviolet cutoff scale r is not too small: 
(rFj )'''»«/' — 1 <C 1]. Of course, in order for the statistical sample to be significant, it is necessary that 
the range of scales r <l <V^ be large enough to accomodate all relevant interactions: (a — 1)T^ /r>l, 
and that enough wavepackets be present in the volume Vi: i.e. { ,J^i\ YVi 3> 1. With a degree of nonlocality 

(a — 1)^^ of at most ten (see [19] and results in the next section), and — in the range of a few percents 
[2], these conditions can be assumed to be satisfied, and all statistical averages, indicated by < >, will be 
understood here to be carried on in space, over a single large but finite volume Vi. 

The closure that is going to be used in the next sections is of the Quasi-Lagrangean type [21], in which 
fluid structures at scale / are studied in a reference frame moving with a speed given by the average fluid 
velocity in a ball of radius //A, A < 1; the free parameter A is then adjusted to lead, in a mean field theory, to 
values of the Kolmogorov constant in agreement with experiments. This mean field theory is obtained using 
still a global Fourier basis. Notice that here l/X is not the wavepacket size, and in this respect, the present 
approach differs from that of Nakano [17]. Of course our approximation is completely uncontrolled, in that 
sweep effects from scales between / and l/\ are still present, while part of the strain from scales larger than 
l/\ is lost. However, it is still an improvement upon using infrared cutoffs in the expression for the energy 
transfer, especially in view of the fact that this last operation would not preserve the nonlocal character of 
the nonlinear interaction and would lead to transfer proflles in disagreement with [19]. 



The fluctuating dynamics is obtained studying the energy equation for the wavepackets [which, after 
angular integration, becomes the energy equation for the shells kn < k < kn+i]- This is obtained by 
substituting the statistical averages which would lead to a mean field closure, with partial spatial averages 
over wavepackets volumes. This leads to an energy equation in the same form (in Fourier space) as the 
original mean field one, plus a noise term that is essentially: < «*"'(«<"' V)w*°' >jj (r,t). Here w'°' is the 
gaussian random field, which gives the lowest order approximation for the fluid velocity. The end result is 
an equation for the energy content of shells at a certain position in space, which is in a form very close to 
the stochastic chains studied by Eggers in [22-23]. 

This paper is organized as follows. The closure technique is going to be described in section II and 
equations for the energy balance in laboratory as well in Lagrangean frame are derived in the mean field 
approximation. In section III, this closure is applied to the analysis of the energy transfer fluctuations. In 
section IV, an energy balance equation in terms of shells is derived and its solution is used to obtain the 
intermittency corrections to scaling. Section V is devoted to discussion of the results and to concluding 
remarks. Most of the calculation details have been confined to the appendices. 

II. Mean energy transfer in Lagrangean reference frame 

A. Closure outline 

In the inertial range of fully developed 3D turbulence, the dynamics obeys the Euler equation: 

^t^r{r, i) + (v • V)v(r, t) + VP(r, t) = 0, (1) 

while the presence of dissipation is accounted for by the energy flux towards small scales; the pressure P is 
calculated through the incompressibility condition V • v = 0. In a Lagrangean reference frame, Eqn. (1) can 
be rewritten in the equivalent form: 

Av(zt +r,t) = ((v(zt, i) - v(zt + r, t)) • V)v(zt +r,t)- VP(zt +r,t) = 0, (2) 

where Dt — dt + v(zt,t) • V is the material derivative along the trajectory zt = zt(ro,to) of a Lagrangean 
tracer that at time to was at Tq. 

The first assumption of the model is that the velocity field can be taken to lowest order to obey gaussian 
statistics, with correlation given, in a Lagrangean frame (with corrections to the exponential decay in t for 
small i), by the expression: 



- 27r2C'Ko;e'/'P"^(k)A:-"/3 exp(-ryfet). (3) 

CaT,(r,i;ri,ti)EE<<'(zt + r,i)w^'"(zt, +ri,ti)>, 

where t > 0, CkoI is the Kolmogorov constant, e the mean energy dissipation, P"'''(k) = 5"^ — '' J^ the 
transverse projector and rjk ~ pl^'^k"^'^ the eddy turnover frequency at scale fc^^, with p a dimensionless 
constant. The statistical average < > taken in Eqn. (3) can be understood as an average over the initial 
position To, which is taken in a large but finite volume, as discussed in the introduction. 

The first nongaussian correction v*^' is obtained from (2); this is integrated, keeping second order terms, 
in the Direct Interaction Approximation (DIA) [24-26], rather than following [16]. The result contains then 
a green function: G{t, r|ii, ri), which gives the effect of a source, which at time ti was in z^^ -l-ri, on a point 
•Lt + '^ ^^ time t: 

f<:'(ro+r,0)- f dt fd^sGPjO,v\t,r + s) 

J — OO J 

X [«^"'(zt + r + s,f)-t;^">(zt,t)]aX'(z*+r + s,t). (4) 
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The time decay of the (retarded) green function introduced above is assumed to be the same as that of the 
correlation function [26] and the foUowing approximation is adopted: 

G(t, r|ii, n) ~ G{t, r - njii, 0) = G(r - n, f - h), (5) 

which means that the initial and final position in G arc shifted together until the first lies on the Lagrangean 
trajectory zt(ro,to)- Eqns. (5) allows a great simplification in taking Fourier transforms, in that all correla- 
tions become diagonal: < VpVq >cx S{p — q). In particular, it becomes possible to write time correlations in 
the form: C^^{t) ^ G;^,fc(i)Cr(0), where: 

G^^(t) = / dVe^'^-G^'^Cr, t) = P"'^(k) cxp{~7jk\t\). (6) 

The meaning of this approximation is that the divergence of Lagrangean trajectories is disregarded. In 
particular, the statistical average over the initial position Tq coincides identically with spatial average at the 
given time. We have also the result that, in this approximation, the advanced green function coincides with 
the "transpose", with respect to the space slots, of the retarded one; in the following sense: for ti > t2, 

<'(zt, +ri,ti)w;'«(zt2 +i'2,t2) >= /"d3sG^(ii,r2+s|i2,r2)G„,(r2 + s-ri,0) 

= J d3sG^(ti, ri|t2, ri + s)G,^(ri + s - r2, 0). (7) 

Further discussion of these points is contained in Appendix A. 

The final assumption is necessary in order to separate sweeping from straining scales, and is that, inside 
averages, one has to carry on the following substitution: 

< ...[v(z + r) - v(z)] • V... >fe^< ...[v(z + r) - w(A, fc)v(z)] • V... >k, (8) 

where w{X, k) is a smoothing operator acting on v(z) by filtering out Fourier modes above Xk. The hypothesis 
here is that, although all velocity components are integrated along a single Lagrangean trajectory Zt(ro,to)i 
when these components are large scale ones, small scale details of zt will not contribute in averages. 

B. Mean field analysis 

In order to study the fluctuation dynamics in the energy transfer T, we will have to work in a Lagrangean 
frame. However, in order to fix the free parameter A introduced above, it is necessary to match theoretical 
predictions with experimental data of the Kolmogorov constant that are taken in a fixed laboratory frame. 
For this reason, the first step in the analysis is the derivation of closure equations for eulerian correlations. 

The calculations to obtain the energy equation in the laboratory frame are standard [26] : by multiplying 
Eqn. (1) by v(ro,0) and taking the average, the pressure term drops off because of incompressibility, and 
one is left with: 

dtC'''>ir,t)\t=o = -< i;"(0,0K(r, 0)9^1;^ (r,0) >, (9) 

which of course is equal to zero at steady state. Actually, this is the equation for the Eulerian 2-time 
correlation at zero time separation; the energy equation is obtained from (9) by multiplying its RHS (right 
hand side) by 2. Expansion to first order in w*^' and use of Eqns. (3-4) and (7), leads to the same expression for 
the energy equation as in DIA (both Eulerian and Lagrangean, [24-26]), and in "Quasi Normal-Markovian" 
closures [20]: 

^^ = (J) V(fc) = ^/ dpdqkpqe,p,b,p,C,{Cp - Gfe), (10) 

where: C'^"' {t) — P"'''(k)Gfe(i); Ck = Gfe(G); A is the domain defined by the triangle inequalities: p > and 
\k — p\ < q < k + p] Okpq = {rjk + rjp + Vq)~^ is the relaxation time; bkpq = {p/k){xy + z^) is the geometric 
factor, in which x, y and z are cosines of the angles opposite respectively to fc, p and g in a triangle with sides 
kpq. The terms associated with integrating along z arc uniform in space, they do not couple with the others 



and arc shown in Appendix A not to contribute to the final result: in this model, the choice of integrating 
back in time along a Lagrangean path is felt only in the eddy turnover time Okpq- 

The green function G is a Lagrangean object and is obtained by multiplying this time Eqn. (2) by 
v(ro,0) and then taking averages. The resulting triplet term is in the form: < f"(ro, 0)[w^(zt + r,t) — 
ti''(zt, i)]9/3W^(zi + r,t) >; now, v^(zt,t), that is the term coming from the shift to a Lagrangean frame, 
contributes to the final expression and is responsible for the cancellation of the sweep terms. Substituting 
Eqns. (3-4) and (7) in the new triplet, and expanding again to first order in w<^', we obtain at steady state 
the following equation: 



DCkjt) _ 1 
Dt 47r2 



dr dpdq[bkpqGp{t - T)Ck{T)Cq{t - 



+^p(A, k)b'klqGp(t - T)Ck(t - T)Cq(T) - wq{\, k)b'-^l^Gu{t - T)Cp(T)Cq(t - r)] (11) 

where: G^^ = P"^Gfe, b^'^^ = {p/2k){xy{l - 2z^) - y^z) and b'^^^ = {p/2k){xy + z(l + z^ - y^j) are new 
geometric terms, and ■Wq{X,k) gives the effect in fc-space of the cutoff operator w{X,k) The term in Eqn. 
(11), which cancels the divergence of the integral for q ^> 0, is the one in &S„- Integrating from i = to 
t = CO and using Eqns. (3) and (6), we get the result: 



r,2 



-^— = - I dpdqpq^^/^ 



CkoI '2 J a 



.b 



'Wr,0,'^ Wa 



kpq . P kpq •? kpq 



p2/3^q2/3 (l+p2/3)q2/3 (1 + ^2/3)^3 



(12) 



which gives a first equation connecting CkoI and p with A. A second equation connecting CkoI and p, given 
an energy balance equation in the form of (10), was derived by Kraichnan [27]: 

p/GL, ^ 0.19. (13) 

The constant CkoI considered in this section is a well defined quantity, provided the average volume Vi and 
the range of scales k are not too large; in this sense we have locally: CkiVi) oc CkoI^i k~^^/^ ^ even though 
for Vi ^> oo anomalous corrections become important, and CkoI ceases to have a clear meaning. Imposing 
that the Kolmogorov constant matches the experimentally observed value CkoI ~ 1-5 (we are assuming here 
that that the finite size effects and intermittency corrections which affect the experimental CkoI are indeed 
small), and using a gaussian profile for the cutoff w: Wp{X,k) = exp(— (p/Afc)^), Eqns. (12) and (13) set: 
A = 0.9. In the following we will fix therefore: 

Wp{X, k) -^ Wp{k) ^ exp(-1.23(p/fc)2). (14) 

The same steps leading from Eqns. (1) to (10) can be repeated starting from Eqn. (2). The result is the 
following equation for velocity correlations in a Lagrangean frame: 



DCk{t) 



Dt 



t=o 



{^'TiX',k) 



1 

47r 



dpdqkpq0kpq[bkpq + Wq{X' , k)b^k^^]Cq{Cp - Ck), (15) 



where: 6jf' ~ (p/fc)(l + xy + z{z^ — ^ ^^ )) [28]. Again, the corresponding energy equation is obtained 
by multiplying the RHS of (15) by 2 and understanding the t on the LHS (left hand side) not as a time 
separation in a 2-time correlation but as the t-dependence in a non stationary 1-time correlation. Notice that 
the parameter A' entering Eqn. (15) does not have to coincide with A — 0.9, which is fixed and separates 
between sweeping and straining scales. This allows to calculate the Kolmogorov constant and the parameter 
p in different reference frames. The result is shown in Fig. 1. and suggests that, although differences 
are expected between quantities measured in laboratory and Lagrangean frame (due to correlation among 
reference frame velocity and quantity to be averaged), their orders of magnitude should be the same. This is 
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particularly important because the intermittency estimates that are going to be derived in the next sections 
will have to be calculated in a Lagrangean, not a laboratory frame. 

The structure of the energy transfer is studied, following Domaradzki and Rogallo, [19] by decomposing 
r(A, k) in its contribution from different scales: 

T(A,fc) ^ f dpT{X,k,p). (16) 

It appears that the present closure is able to maintain the features of large scale straining observed in [19] 
also in a Lagrangean frame, as it should be expected. It is clear instead, from Fig. 2., that simpler closures 
based on a Navier-Stokes nonlinearity in fc-space, amputated of the large scale convection contributions: 

(fc„P^^(k) + k^Pp^{k))v^v^^ ^ (fc„(l - w,)Pfj-,{k) + fc^(l - Wp)Pfj^{k)))v^vl (17) 

would lead to transfer profiles almost without any exchange of energy between nearby scales. 

III. Energy transfer fluctuations 

Although the analysis in section II disregarded fluctuations, all averages were implicitly dependent on 
the position (through e), at the scale of the volumes Vi. In this section, we calculate the same averages 
over balls of radius R and the variations of the result from ball to ball are used to study the fluctuations of 

— 1/3 . 

the energy transfer. If the ratio RV^ is not too large, these fluctuations are going to be small and the 
analysis can be carried on in perturbation theory. 

In principle, it would be nice to carry on the calculations in the laboratory frame, where experimental 
data are available. However, the statistics of the velocity field is ill defined there, due to the effect of sweep. 
The balls are then imagined to move rigidly along Lagrangean trajectories passing through their centers 
at time t^, while, thanks to the simplifying assumptions of Eqns. (6) and (7), the average over the initial 
position is substituted by a spatial average inside the balls. As discussed in the introduction, we associate 
to each scale /, a certain radius R{1) ~ l/{a — 1); this automatically induces a basis of wavepackets of 
width Afc ~ R~^ . In particular, we obtain a partition of Fourier space in shells S„ of radii fc„ = a"fco, 
each containing ^ (fc„/Afc)^ wavepackets, associated with the different orientations of the wavevector k [17]. 
Clearly, a further degeneracy in the wavepackets is produced by their different location in real space. 

For R/l large, the main interactions occur among overlapping balls and are associated with local transfer 
of energy towards small scales. We consider therefore a sequence of nested balls and study the transfer of 
energy among them. Notice that a derivation of deterministic equations for a shell model of the type 
considered here, would require some dynamical analog of the statistical assumptions of Eqns. (5-7); at this 
point it is more natural to follow the route of statistical closure to the end. 

Let us indicate with < >„ the spatial average over B^: the ball of radius i?„i associated with the ?Ti-th 
shell; once we have fixed the origin of the axis of the moving frame on z^, we can write < >,„ in terms of a 
kernel VF(r, m): 



In this way, the energy density in Bm reads: 



//■ d fc 



E{m)^ j dkEk{m)- Ek{m) ^ 2TTk^ J -^Wp{m)w^+^ ■ v^_^, (18) 

To fix our ideas we shall consider gaussian wavepackets: 

Wk{m) = exp(-(fc/Afc„)2), Afc„ = c^{a - l)fc,„, (19) 

with c^ relating the wavepacket thickness Afc„i and the shell spacing km+i — fcm = (^ ^ l)fcm- Notice that 
Cw can be treated as a free parameter, in that, for a fixed choice of shell (i.e. for a given value of a), we can 
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still consider arbitrarily thin wavepackets, or in the limit, even global Fourier modes. Of course, it is the 
thickest wavepacket for a given a (that is the maximum value of c^,), which will be able to catch most of the 
fluctuation dynamics. 

If knRm is large, Wp(ra) ex R^5{'p) and we can approximate: 

Ek{m) ~ (1 + 0„(to))4"' = (1 + c^,,{m))CKoil^'^k-'>'\ (20) 

with £'^"' = Cfefc^/27r^ the spectral energy density in Vi, n — n(k) the shell of the wavevector k, and 
(f'ni'm) = 4'n{m,t) fluctuating and small. 

The energy density E{m) can be expressed also as a sum of contributions from different shells: 

E{m) = Yl En{m)- En{m) = f dkEk{m). (21) 

The term En{m) ~< E(k^^, a, Zt, t) >m can be seen as the istantaneous total energy of wavepackets in S„, 
lying in the volume B^ at z^; however, if B„i becomes smaller than the wavepackets, i.e. (for c^ maximal) 
n > m, the only average taking place will be over wavevector orientations and will be independent of to; 
hence we get: (f'ni'm < n) = 0n(n). 

With these definitions, the energy equation for the shells reads: 

DtEnim) = T,,{m) + /„(m). (22) 

In the equation above, the energy transfer into S„, averaged over B„i, has two components. The first: 

r„(TO) = I dkT'ik), (23) 

J n 

is a relaxation term giving the response of the system to fluctuations in the energy content of the various 
shells. For small (f), the integrand T'(fc) coincides with the transfer term T{X,k) of Eqn. (15), with the 
substitution: Cfc ^- (1 + (j)n{'m))Ck- For c^ — > c^^^ , the dfc integral in Eqn. (23) receives discreteness 
corrections from the shell and the wavepacket thickness being comparable; notice however that the J. dpdq 
integral contained in T'(fc) [see Eqns. (10) and (15)] remains unaffected, when a — 1 is small. 

The second term contains the fluctuations of the transfer, around the average. To lowest order in 
the expansion around gaussian statistics it has itself two components. The flrst comes from the pressure: 
< V • VP >m and results in a surface integral over the boundary of B^. The second: < v<"> • (v<"> • V)v<"> >„ 
did not contribute in the mean field analysis of the previous section, because of the gaussianity of v*"'; it 
plays a role here however, by acting as a source of fluctuations in Eqn. (22). Now, we are considering a 
situation in which kR is large; hence, the pressure contribution can be neglected and we are left with: 

/„(m) = - < h{n)[w^°\zt,t) ■ ((v<">(z, +r,i) - v<'"(z,,i)) • V)v<">(z, + r, i)] >„, (24) 

where h{ri) is a band pass filter for the modes in S„: h{ri)f{r) = J d^r'/i(n, r')/(r + r') = J ,^3j fefc("-)/ke''^''', 
with hk{n) = H{k — kn) — H{k — fc„+i); H{x) is the Heaviside step function. 

Treating fn{in) as an external (nongaussian) noise, Eqn. (22) becomes a stochastic differential equation 
for the (j>n{m)'s. Following Eggers [22-23], we express fn{iTi) as the difference between the fluctuations in 
the energy flux across kn and kn+i- 

fn{m) ^ gn(m) - gn+i(m), (25) 

where g„ and gn+i are associated one with each of the step functions entering the deflnition of hk{n). 

A. Fluctuation source 

The statistics of 5„(to.) can be calculate explicitly from Eqns. (3) (24) and (25). Here we shall content 
ourselves with the analysis of the 2-point correlations. Already, this calculation requires the evaluation of 
some fifteen contractions of the product: 

w"(zi)[«'^(zi + n) - «''(zi)]a;3«"(zi + ri)w^(z2)[«"(z2 + ra) - v%Z2)]daV^ {^2 + ra); (26) 



these contractions are listed in tabic Bl. in Appendix B. Of these only six contribute; one example is 
contraction (B1.9): 

< .g„i(TOi,ti)g„2(TO2,t2) ><"> = 



d^rid^r2i/(ri,ni)if(r2,n2) / d^Zid^Z2VK(zi, mi)W^(z2, m2)C"''(zi - Z2, f) 
xa,[C^T(zi + n - Z2 - r2, t) - 7«C^T(zi - Z2 - r2, t)] 



X 



a^[C""(zi + n - Z2 - r2, t) - 7«C""(zi + n - Z2, t)]; (27) 



notice the cutoff operators li signaling a term coining from working in a Lagrangean reference frame. In 
terms of Fourier components, we get: 

< 9ni{mi,ti)gn^{m2,t2) >'"' = 

d^fc A^p A^q 



rCkCpCqCxp{-{r]k +rjp + riq)\t\) 



(27r)3 (27r)3 {2Tif 

xWp+q-k(™i)VKp+q_k(m2)pfcyz(j/ + a:;z) 

X [(2 - Wp{k))H{k - fc„J - Wp{k)H{q - fc„J] 

X[(2 - Wq{k))H{k - fc„J - Wq{k)H{p - fc„J]. (28) 

(The wavevectors entering the two last lines of Eqn. (28) can be tracked back to Eqn. (26) and (B1.9): the 
p and k in 'Wp{k) come from v^ and w^, while the k and q in H{k — fc„j) and H{q — km) come from w^ and 
v"). It has already been mentioned that, for large kR, the averaging kernel W is proportional in fc-space to 
a Dirac delta. Also the product W {mi)W (1712) is a Dirac delta: 

W^{m,)W^{m2) o, (cUa-iyf(^^^^y^^S(k), (29) 

where use has been made of Eqn. (19). Eqn. (29) allows simplification of Eqn. (28) by means of the bipolar 
integral formula [26]: / d^pd^qS{p + q — k) — -^^^ J^ dpdq. Repeating the calculations leading to (28) with 
all the other contractions we obtain the result: 

< gni{m,ti)gn2{m,t2) >= — —[CkoICw (a- l)km] 

22 
/•oo /• 

X / dk dpdq{kpq)-^/^cxp(-{'i]k+Vp + Vq)\t\) 

Jo J A 

x[HiBi{kpq) + H2B2{kpq) + HsBsikpq)], (30) 

The terms Bi i = 1,2,3 are geometrical factors similar to the bkpq terms entering the expression for the 
transfer function. The factors Hi = Hi{ni^2', k,p, q) are expressed in terms of step functions and restrict the 
integrals to the appropriate domains for the calculation of Lagrangean frame energy fluxes; they have the 
same origin as the two last lines of Eqn. (28). The exact form of both the Bi and Hi functions is given in 
Eqns. (B4-7) of Appendix B. 

For i = 0, Eqn. (30) can be reduced to double integrals using the standard change of variables [26]: 
k — kn/u, p — knv/u, q — knw/u and exploiting the similarity of the integrand with respect to u. Similarly, 
the following expression for the noise correlation time: 



77„(m) =< fnimy > / dt< /„(to, t)/„(m, 0) >, 

Jo 

can be reduced, after explicit calculation of the time integral, to a double integral. Numerical evaluation 
gives then the results: 

< gn{mf >~ 0.15 {CkoI c^ {a - l))-''a'("^")e2. (3^) 



and: 

^„(m) - pei/3fc2/3 p~0.6, (32) 

where l{n,m) = min(0,TO — n) for c„, maximal. The decay of the correlation with respect to scale is 
shown in Fig. 3. Notice that p ~ 2px=a.9 (see Fig. 1.), which is what one would get identifying naively: 
< g{Q)g{t) >oi< v{0)v{t) >^. Similarly, the decay of < gnQn' > with respect to (n — n') In a is approximately 
the same as that of the auto-correlation for |T(fc,p)|: A{q) — J dp\T{k,p)T{k + q,p + q)\, which is due to 
the fact that the energy flux at scale k is roughly equal to J dp\T{k,p)\. 

B. Relaxation term 

The shell energy equation (22) is in the form of a nonlinear stochastic differential equation. Since (j) is 

3 i 

small, we can linearize Tn(m) (including the frequencies rjk ~ ^^E^ entering the term Okpq), with a result 
that is similar to the stochastic model of Eggers [22] . One difference is the presence here of energy transfer 
between non adjacent shells. To analyze this issue, we write the transfer Tn{m) as a sum of contributions in 
which the wavevectors k p and q are respectively in S„, Sr and S^: 






r,s 



where, from Eqns. (15) and (20): 

Cnrs = j '^^ j dpdqOkpq akpqk^ipqY^^'^ - -P^^^Okpqicikpqk^ipqy^^^ 



Ar 



-bkpqp'^ikq)-''/^ - bkqpq^{kp)-^"'^ + bkpqq^ikp)-''"-^, (34) 

and: 



dnrs = I <ik dpdq Okpq 



^|fc'/'4p,(afepgfc'(p<z)"''/' 



~bkpqp\kq)-'''^ - bkqpq\kp)-'''^^ + hpqp'ikq)-'''^ + 5fe,pg2(^p)-5/3J _ (35) 

The quantities appearing on the RHS of Eqns. (33-34) arc defined as follows: A^s is the restriction of A to 
p e S^; g e Ss; bkpq = bkpq -h Wq(k) b'-j^^^; hkpq ^ bkpq + bkqp and hpq == (fc^/^ +p2/3 _^ g2/3)-i^ rj,^^ j^jjg ^^ 
Eqns. (33-34) can be reduced to double integrals using the same method of Eqn. (30), and are evaluated 
numerically. We can rearrange the sum in Eqn. (33) in the following form: 

00 

r^s r r— 1 

where A is the finite difference operator acting as follows: 

A2'-+V„ = ^(A2''0„+i - A2'- </,„_!); 

A2>„ ^ A2(--i)0„+i + A2(--i)0„_i - 2A2('-i)0„, (37) 

so that the coefficients A^ (m) are expressed in terms of the a^ {m) 's through the relation: 

A2r+l — an+2r-l — 0,n-2r+l', ^2r — ;: (oo) 

The fact that only differences between </>„ enter the expression above is due to the fact that we are expanding 
around a Kolmogorov spectrum, for which T„ = 0. For the same reason, the coefficients Cnrs and dnrs are 
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invariant under the transformation: {n, r, s} ^ {n + j, r + j,s + j}, which explains the fact that the finite 
difference coefficients A^ do not depend on the sheU index n. 

In Fig. 4 these coefficients are plotted in terms of the shell constant a. Notice that finite difference of 
order higher than 2 appear to be negligible for most choices of a. We obtain therefore the basic result that 
the shell dynamics obeys, for most values of a, a (discrete) heat equation forced by a random noise, with an 
advection term proportional to A. 

The physical picture corresponding to this result is not new [29]. Fluctuations over finite volumes 
Bm moving with the flow, in the energy transfer to eddies at scale fc„, generate fluctuations in the energy 
content of these volumes and scales. The eddies are stretched by the turbulent flow, so that their energy 
(together with its fluctuation 0„) is transferred towards smaller scales; the term responsible for this effect 
is the advection ^i A0„. At the same time, the randomness of the turbulent flow causes some eddies to be 
stretched more and some less, resulting in the end in a diffusion of energy over different scales. 

IV. Anomalous scaling estimates 

The fluctuations in the energy of wavepackets at different space locations is what is responsible for 
intermittency in this model. Instead of studying the scaling behaviors of the structure functions S'(Z, q) —< 

^1 



V? >, we focus on the modified structure functions, defined through: 



Sik-\a,c^,2q) ^« En >l>^ Ei°> <(l + (/.„(n))«>. (39) 

Identifying naively the space separation I in S{1, q) with the wavelength k~^ in Eqn. (39), the two definitions 
of structure function coincide, for an appropriate choice of a, and for c^ = c^^"*"^. Given a large enough 
averaging volume Vi, the intermittency correction is small and we can linearize in n; in this way, using also: 
kn — koa"', we have: 

(:^^-qC,c--^^<{l + Mn)y>- (40) 

in a an 

If we confine ourselves to the lowest order moments, large fluctuations of (j) do not contribute too much, so 

that we can expand: (1 + (/)„(n))'' ~ 1 + "^^"^^ ' < (/)„(n)^ >. Imposing the Kolmogorov relation for the third 

order structure function: 5'(/,a,3) ex /, fixes the value for 6(^2, leading to the lognormal statistics result [6]: 

we see then that the presence of anomalous scaling is associated with secular behavior of the fluctuations 
(f>„{n). 

Unfortunately, the energy balance for the wavepackets, Eqn. (22) is not an equation for (f>n{n), but one 
for 0„(to) for fixed m. We can rewrite Eqn. (22) in a more explicit form, by approximating the relaxation 
term Tn(rn) with the advection-diffusion operator introduced through Eqns. (35-37), and by rewriting the 
noise term /„ as a difference of energy flux fluctuations at different scales, as from Eqn. (25). After adequate 
rescaling, we are left with the equation: 

[exp{^-fn)di - DA^ + FA)]0„(TO,t) = A(F„(m)i/2e), (42) 

where A^„ = ^„+i - ^„ and: 

- ,__l/3,2/3 2 CkoI . 

V = _.^-°' . 4,; F . 0.83 -(^-('^ - ^»^^-' «""'"^ ■ 



2 



pp(l-a-2/3) ^' - ■ (p(l-a-2/3)) 
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Here the Kolmogorov and time scale constants CkoI and p are the ones measured in the moving frame: 
CkoI — 0.8 and p ~ 0.3 (it appears however that all final results do not change by more than 10% by 
exchanging these values with their Eulcrian counterparts: CkoI ~ 1.5 and p ~ 0.43). The dependence of the 
noise correlation < ^^ > on n — n' in the formula above is a fit of the result of numerical integration shown 
in Fig. 3. 

Clearly, Eqn. (42) does not lead to intermittent behaviors; the factor a'^™'"' — a^™^"' in the noise 
amplitude, which goes to zero at small scales prevents it. Passing from Eqn. (42) to an equation for (t>n{n) 
requires the introduction of corrections due to the fact that now, wavepackets associated with different scales 
do not overlap exactly in real space. There are two such contributions: 

{-DA'' + yA)(/.„(m)U=„ - (-Z5A2 + yA)(/.„(n) 

= {-D + l//2)(0„+i(n + 1) - ci>n+i{n)) (44a) 

and 

k(F^/'{n, m)C„(m))|„=„ - k{F^/'{n, n)Un)) = 

F^/'{n,n){Ca+i{n + l)-^n+i{n)). (446) 

The shell equation for (/>„ = 0„(ri,i) takes then the form (to simplify notations, the hat on the rescaled time 
t will be dropped in the following): 

[exp(-7n)9t - DA' + VA)]cf>^ = F^/\M + SO, (45) 

with 60 which is equal to the sum of the RHS's of Eqns. (43a-b), providing a new fluctuation source beside 
the original term A^. 

At this point we are in the condition to identify the terms in the shell energy equation which are 
responsible for the generation of intermittency. For the sake of simplicity, and to make contact with the 
model of Eggers [22], let us adopt for a moment the approximation: O — e"'''"/^^ with < ^(t)^(i') >= S{t—t'), 
and similar equation for §0 The RHS of Eqn. (45) is then proportional to: 

AC -7^/2 + -56 (46) 

We see then that Eqn. (45), at statistical equilibrium, is very similar to a random walk equation in which 
the role of the time is played by n/V and that of the random kicking by —7^/2 + 60 This noise term 
continuously pumps into the system fluctuations, which are dissipated at very large n, by viscosity. It is the 
random walk character of the process that leads to the linear growth of < 0^ >, with respect to n, already 
observed in [22]. The two terms providing the source of intermittency have different physical origin. The 
term SO which comes from Eqns. (44a, b), is due to the competition in the energy transfer between eddies 
at different locations, characteristic of the Random Beta Model [8]. The term 7^/2 instead, comes from the 
mismatch in the characteristic time scales of the energy transfer between different shells and was responsible 
for the production of intermittency in the stochastic model of Eggers [22-23] . 

The first term in Eqn. (46), which is the derivative of a random noise, produces the gaussian part of 
the fluctuations in the energy content of the wavepackets. Notice that this same quantity can be calculated 
directly from finite volume averages of |v<">|^, given the expression for the correlation given by Eqn. (3), 
together with Eqn. (18). This fact will be used to provide a check on the goodness of the approximations used 
to arrive to equation (45). Notice finally how in this approach, the smallness of the anomalous corrections is 
associated with the smallness of the parameter a— 1 and with the fact that the amplitude of the intermittency 
source term is second order in this quantity, with respect to the source term of the gaussian fiuctuations. 

A. Solution of the shell energy equation 

It is possible to solve Eqn. (45) either following [22], by diagonalizing the LHS of (45), considered as a 
matrix equation, or using the multiplier technique adopted in [23]. Here, the smallness of the parameter a— 1 
allows to consider the continuous limit of Eqn. (45) and to use a multiple scale expansion in which, to lowest 
order, the dependence on n produced by the exp(7n) terms is neglected on the scale of the fluctuations. This 
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allows a solution of the problem in terms of green functions, in which no approximation on the form of the 
noise correlation is required. The green function for Eqn. (45) is: 



exp 



{n—n^ — V e'^^t) . ^n 



The first quantity that we are going to calculate is the gaussian part of the fluctuations: 

<^%>=Fj^^\gUn)\'k' <\^k.\'> 
^F^[ dk (l + i?P)exp(-M4|fc|) 

where gfcaj (jT-) — {—i{uje^'^"--Vk) + Dk'^)^^ is the Fourier transform with respect to n' andt of g{n,n' ,t). As 
mentioned before, we can repeat the calculation directly from the statistics of the velocity field v'°'. Consid- 
ering the limit of Cw small, corresponding to wavepackets thinner than the shell, nondiagonal contributions 
in Eqn. (18) can be neglected and we have: 



=.0.032"-°'""°-.'°-"'°. (49) 

(1 — a 3 )2 

The two expressions for < (J>q > given by Eqns. (48) and (49) are plotted in Fig. 5. against a; the best 
agreement, though still rather rough, is obtained for the range 1.2 < a < 1.3, which is consistent with what 
would be expected by looking at the energy transfer profile (Fig. 2.) and at the one for the correlation of 
the energy flux fluctuations (Fig. 3.). In the same way it is possible to calculate the correlation time for (J>g- 

- V-T / dt < 0G(O)0G(t) > 



<rG> h 
-iF f dfc duj fc2 



<(pn> 



/dfc duj fc , , , ,9 , ,9 , , 



The prediction from direct calculation, obtained from the generalization of Eqn. (49) to 2-time correlations, 
is Tc(fco) — Ij corresponding, in non rescaled units, to Tc{k)~^ = 2r]k; here we find that 0.9 < rc(fco) < 1 in 
the whole range 1.1 < a < 2. 

We turn next to the calculation of the intermittent part of the fluctuation (f>. Care must be taken now 
due to the divergent nature of correlations at large n, which forbids in particular the use of the Fourier 
representation adopted in Eqn. (48). 

The equation giving the growth of < cf)'^ > at large n is obtained by multiplying Eqn. (45) by and 
taking the average: 

< 0„[exp(-7n)at - DA'' + VA]q^^ >= F^/^ < 0„(Ae„ + S^ > ■ (51) 

We subtract from Eqn. (51) the gaussian part of the fluctuations as given by Eqn. (48); then the flnite 
difference A acts only on the part of the noise variation due to the scaling of the correlation time, which 
is of order a — 1. Next, we notice that to lowest order, the contribution to (5^ coming from Eqn. (43a) is 
obtained by approximating (f) by its gaussian component (/)g. The RHS of Eqn. (51) takes then the form: 

S(a, c„) - F /" dt I dn' [g{n, n' , t) (7^(1 - t") + /?) 
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+.g(n,n' - l,i)/3'i] < e„(0)e„+„'(t) > +/3 < 0gP - y^ 



/>G >, (52) 



where (3 and (3' are ©((a — 1)^) quantities giving respectively the amount of wavepacket volumes not over- 
lapping (Beta model effect), and the correlation between this effect and that of the scale dependence of the 
noise correlation time. Explicit expressions for these coefficients and further simplification of Eqn. (52) are 
presented in Appendix C. 

Next turn to the LHS of Eqn. (51). The time derivative term is equal to zero at steady state, while it 
is shown in Appendix C that the intermittent part of the fluctuations does not contribute to the < (j)A'^(j) > 
term. We are left then with: < (j)A(j) >~ ^A < (jp >, so that we obtain the result, for the kurtosis scaling 
exponent: 

or ^ ~ ^ d<^^> 2S(a,c^) 
in a dn V (a) In a 

The dependence of 2^2 — C4 on a for fixed c^ is shown in Fig. 6.; notice the saturation occurring at a ~ 1.3, 
suggesting that the bulk of intermittency production occurs at scales of the order of three to ten times the 
size of the eddies in exam. From inspection of Eqns. (43), (52), (53) and (C3-6), we see that, for small a— 1: 
2C2 — C4 = ^{(^wi^^ ~ l)"^); so that we obtain a direct connection between the smallness of the intermittency 
corrections and their being proportional to a rather large power of the small parameter a — 1. Unfortunately, 
the intermittency exponents defined here depend sensitively on the parameter c^,; due to the difficulty in 
determining c'^^^ with precisions, it is therefore problematic to make accurate predictions on the scaling of 
S{1, q) starting from S{1, a, Cw, q)- 

A value 2^2 — C4 — 0.016, is obtained setting c„, — 2, corresponding to the "reasonable" condition that 
the noise granularity be equal to the shell thickness (a — l)fc„, i.e., from Eqns. (19) and (29): Wp^^_y^ ~ 

exp(— ^J^^'^"^^ y ). However choices as reasonable as the one just considered, lead to results differing from 
one another by up to an order of magnitude, so that the above formula should not be taken too seriously. 

V. Summary and conclusions 

We have described a model for the production of intermittency in the incrtial range of three dimensional 
turbulence, based on statistical closure of the Navier Stokes equation. A connection between Navier Stokes 
dynamics and phenomenological models like the Random Beta Model [8] , and the stochastic chains considered 
in [22-23], is in this way established. Although this connection may be rather tenuous, because of the 
assumptions adopted in deriving the closure, it is still pleasing that the results in this paper are obtained as 
lowest order corrections to a mean field approximation, which by itself would produce Kolmogorov scaling. It 
should also be mentioned that the present approach is able to produce dynamically, values of non-intermittent 
part of the energy fluctuations in agreement with the prediction from quasigaussian statistics of the velocity 
field; this is a bonus, which provides an indirect check on the goodness of the closure technique adopted. 

The use of wavepackets rather than Fourier modes, is the reason why a perturbative treatment of 
intermittency has been possible here. There has been indeed an intriguing aspect in this subject, namely, 
the difficulty in associating, to the smallness of the anomalous corrections, a small parameter in which to carry 
on perturbation theory. One of the results of the present model is the identification of this parameter with the 
physical quantity Afc/fc, i.e. the ratio between eddy size and the scale of the characteristic flow straining the 
eddy. Scaling corrections appear to be of third order in this quantity, with equally important contributions 
from Random Beta Model kind of effects [8], and from the mechanism of intermittency production of the 
model studied by Eggers [22-23]. 

The main result of this paper justifies a posteriori the use of wavepackets in the problem: the saturation 
in the a-depcndence of the generalized structure functions S{1, a, c^, q), for a > 1.3, which is consistent with 
a separation of scales between production of intermittency and energy transfer. This is a definite prediction 
of the model, which, together with predictions on the actual magnitude of scaling corrections, could be tested 
by direct analysis of experimental data in terms of the generalized structure functions S'(/, a, c^, g). This 
would extend similar studies, carried on by Menevau [30] using wavelet analysis. 

We have not been able yet to establish a quantitatively accurate relation between the two structure 
function definitions: S{l,a,Cw,q) and S{l,q), the reason being the sensitive dependence of S{l,a,c-u],q) on 
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Cm> = / „i-|t. i the ratio of wavepacket to shell thickness in fc-space. Preliminary estimates hint towards a 
value for the scaling corrections smaller by a factor of the order of five than the experimental one. In the 
lognormal approximation: /_g.. ^ —.002, which should be compared with the lognormal fit of experimental 

values: , _^-. ~ —.01 [2]. Due to uncertainties in the present analysis, it is not possible yet to state whether 
incrtial range processes are as important for intermittency production as finite size corrections, or if they are 
themselves just corrections to dominant finite size effects. We reiterate however that such an answer could 
be obtained analyzing experimental data in terms of S{1, a, Cw, q) instead of S{1, q). 

Turning our attention to more formal issues, it is interesting to look for similarities between the model 
described in this paper and the various phenomenological approaches that have been used to study infinite 
Reynolds number intermittency. The basic ingredient here is the partition of Fourier space in shells of 
exponentially increasing radii. For this reason, our approach has a lot in common with the deterministic 
Shell Models studied, among others, by Yamada and Ohkitani [31] and by Jensen, Paladin and Vulpiani [32]. 
An interesting interpretation of these models, particularly clear in Zimin [33], comes from looking at the 
dynamical variable in each shell as the velocity of nested eddies, all located at the same space position, and 
interacting with one another locally in both scale and space. In our approach, such dramatic restriction of 
phase space does not take place, in that the behavior of the various wavepackets in a given shell is treated in 
an average sense; it is then possible that our model may underestimate the amount of intermittency produced 
in the inertial range. It remains to be seen how important this difference is; in principle one should compare 
the present model with some closure for the Yamada-Ohkitani model or, viceversa, look for a deterministic 
dynamical system, whose closure coincide with the one described in this paper, and then compare with the 
Yamada-Okhitani system. 

The connection with the stochastic model of Eggers [22-23] is clearly simpler. In both cases one ends 
up with the same stochastic chain, the only difference being here the not exact conservation of energy 
transferred between shells. If this effect, due to the not exact overlapping of wavepackets at different scales, 
were neglected, the present model and the one of Eggers would coincide. 

Going back to the possibility of turbulent Navier Stokes dynamics taking place in a very limited region 
of phase space, this constitutes the main conceptual limitation of our model. The possibility of preferential 
transfer of energy between individual wavepackets would correspond to intermittency being associated with 
locally anisotropic fluctuations. It should be mentioned that this looks somewhat unlikely, due to the 
separation of scales described earlier and the consequent large number of eddies in one straining region Sk^^. 
In any case this is a possibility, which could lead, in a suggestive way, to a mechanism for the creation of 
coherent structures. In [34], She derived anomalous exponents, based on assumptions on the shape of the 
most intense and rare dissipation regions, so that intermittency and coherent structures may be connected. 
It is difficult to see how shell models, of either the deterministic or the stochastic variety, could be used to 
study such coherent structures, given the very distant scales involved in their dynamics. In any case, perhaps 
this is not a big problen, since, although vortex tubes have been observed in many numerical simulations [3], 
they do not seem to contribute much to energy dissipation and their relevance to the turbulent dynamics 
has been questioned recently in [35]. 
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Appendix A. Quasi Lagrangean closure 

The basic difficulty in dealing with Quasi Lagrangean closures is that studying the Navier Stokes dy- 
namics in a reference frame moving with a single Lagrangean tracer, neglects the divergence of trajectories 
of different tracers. In a more refined closure scheme, this effect would be included by associating each point 
in a correlation function to a different tracer trajectory, like in the Lagrangean History DIA (LHDIA) of 
Kraichnan [25]. However, also in that theory, abridgements of the closure equations were necessary in the 
end, which were similar in their effect to neglecting the divergence of Lagrangean trajectories. 

We try here to implement this approximation in as much a consistent way as possible. Notice first that 
correlations are defined starting from an initial case in which the initial point lies on the Lagrangean trajectory 
Zt, so that Eqn. (3) basically describes the decorrelation of points at distance r from the Lagrangean 
trajectory z^, with respect to points lying on it. Since decorrelation occurs forward in time, this is the 
optimal choice; choosing Z( as the final point and Zq + r as the initial one, would lead in particular to no 
decorrelation, due to the advection term being zero in Eqn. (2). 

The approximation of Eqn. (5) in which the initial point z^^ + ri of the green function G(t, r\ti, ri) is 
shifted on the Lagrangean trajectory, is justified with the previous choice. We try to make it more appealing 
by giving the next order in the expansion around G{t, r\ti, ri) = G(r — ri, i — ti). We have first: 

/"dVe-*"-<t;"(ro,0)w^(zt(ro,0) + r,i) >= /dVe-*"" /d^s G^''(0,r + s|i,r)C;(r - s) 

= / dV I ^GZir, t)C£,, e^(''-''>, (Al) 

where G^''(r,t) is the Fourier transform of G^''(0,r + s|t, r) with respect to s. Expanding the argument of 
the integral in the last line of Eqn. (Al) in k — ki and r, we get then: 

d-V e-'^'' < w"(ro, 0)wT(zt(ro, 0) + r, t) > 
GriO,t)C^.,-l4:r4^^^\Gr{r,t)CS,] +... (A2) 



r=0 



k V ' ^ k,p 2 dk'f' dk'f' dr^ dr^ 
More in general we would obtain: 

dVe-*"- < v'^izt, (ro, 0) + ri,h))v^{zt,+t{ro,0) + ri+r,h+t)> 
Gr(ri,t)C," -iiiA^rG-(ri+r,OG,"J + ... (A3) 



2 dk'f' dk'f' dr^ dr^ 



r=0 



which can be then be expanded in Taylor series around ri = 0. This shows that the lowest order approxi- 
mation adopted in Eqn. (5) is rather bad and that our closure could not be expected to lead to quantitative 
accurate results. 

We turn next to the derivation of the energy balance equation in the laboratory frame, Eqn. (10). We 
consider just one term: 

-T^y'dVidVadVae-'^'^'-^-P'-^-'i'--^) < v^''''-{ruOWp\r2,0)d('v'^^r3,0) > 

= — i-- /dVidV2dV3e-*('^"'i-P'"=^-i""-^) fd^s f drG£(s,-T) 
(27r)-^ J J J_oo 

X < ^"(ri, 0)v^(r2, 0)a^K(z, + rg + s, r) - wi;"(z„ T)]d,Vp{zr + rg + s, r) > . (A4) 

Splitting the 4-point correlation into 2-point ones, we obtain, in terms of Fourier transforms: 

2r5(k - P - q) / dVzdVa J dc'kidc'k2d^k36ik + p + ka - ]^3)0k^k,k,Ck, C'k, 
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where: 



and 



: [B4{kik2k3)[S{ki - k2 - k3)(5(k + q + 2k2 + k^) - wS{ki - k^jSik + q + ka + kg)] 

+B2{kik2k3)[S{ki - k2 - k3)(5(k + q + 2k2+ kg) - W<5(ki - k2)5{k + q + 2k2)]}, (A5) 

B2iklk2k3) = P"''(fci)P„p(fe)fcf P,3a(fc3)fc2 = fclfc2(l + ^')(^ + Xy), (A6) 

Bi{kik2k3) = fc2"Paa(A:i)P"''(fc2)Pp/3(k3)fcf = -^1^2x2(2; + yz). (A7) 



The various terms in Eqn. (A5) are tracked back to Eqn. (A4) as follows: the terms in P4 come from 
the contraction < v"{v" — wv") > d' G^ < vpd,yVp >; those in P2, from the remaining contraction; the 
terms in w come from the shift to Lagrangean frame. Notice now that, from Eqns. (A6-7): 52(^1^1^3) = 
Bi{kik2ki) = 0. Thus the terms in w in Eqn. (A5) disappear and one is left with the same expression 
that would be obtained from Eulcrian DIA [24], but with the Lagrangean response time 0. Repeating the 
same calculation with the other two choices for v''^\ we see that no terms in w contribute and we obtain the 
standard result of Eqn. (10). 

We pass to the calculation of the 2-time Lagrangean correlation: at steady state the part of the time 
integrals from r < do not contribute and we have: 



dV e-^"^"- < f "(ro, 0)K(zt + r,t)- wv^i^^mpVo^i^t + v,t)> 



= I A^re-'^' / dr / d3s{G^( 



S,f-T 



[C"-(-r - s, r) - ^C"^(0, T)]d„dpCP^{-s, t - r) 

+9^[C(-s, i - r) - wCl{v, t - T)]d„CP'^{-v - s, r) 
u.G^(s, t - t) [[C""(-s, r) - u-C""(0, r)]a,a^C^(r - s, t - r) 
+9^[C(r - s, t - r) - iZ;C^(r, t - t)]9,C"'"(-s, r) 

+diiGc.p [[C"" (-r - s, r) - wC"" (0, r)] 
xa,[C^^(-s, t-r)- C^l(-v - s, t - r)] 



+d„C'^P{-v - s, r)[C'=''"(-s, t-T)- wC^^^ir, t - t) 
-wC'''"(-r-s,t-r)+u;wC'3'"(0,t-r)] |. 



(A8) 



Of the terms on the RHS of Eqn. (A8), there is a group which does not depend on the integration variable 
s; these terms come from integrating along the Lagrangean trajectory and lead, after Fourier transform, to 
triads in which one of the wavevectors is zero; as before they do not contribute to the final result. The 
remaining terms in w come from the w on the LHS of the equation, which express the fact that we are 
dealing with a Lagrangean correlation. These terms remain and are responsible for the cancellation of the 
sweep part of the correlation decay. After Fourier transorm, we obtain the result: 



dr 



d?p A^q 



(5(k-p-q) 



(27r)3 (27r)3 

X B3{kpq)Gp{t - T)Ck{T)C,{t - r) - WpB3{q, -p, k)Gp{t - T)Ck{t - T)Cg{T)] 
+ [Bi{kpq)Gp{t - T)Ck{T)C,{t - r) - WpBi{q,-p,k)Gp{t - T)Ck{t - t)C,{t)] 

+ [Bi{kpq)Gp{t - T)Ck{T)Cq{t - r) - WqBiip, k, -q)Gk{t - T)Cp{T)Cq{t - t)] 
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+B2{kpq)G.p{t - T)Ck{T)Cq{t - r) - WqB2{p, k, -q)Gk{t - T)Cp{T)Cq{t - t)]J , (A9) 

where: 

Biikpq) = p,P-''{q)P^p{k)PP^{p)k(3 = ^kpyziy + xz) (AlO) 

and 

Baikpq) = kpP^Pip)Pp^{q)P"^{k)p^ - kpxy{l - z^). (All) 

(The same notation of Leslie [26] is used here for the functions Bi, i = 1, ...4). The terms on the RHS of 
Eqn. (A9) are ordered as on the RHS of Eqn. (A8), once the terms of the last one, that are equal to zero, 
are eliminated. Using the following relations: 

Bsiq, -p, k) - -Biikpq); B^{q, ~p, k) = -B^{kpq) 

Bi{p,k,-q) = Bi{kpq); B2{p,k,-q) = B2{kpq), (Al2) 

and the definitions: 

^f^P" == i E Mkpq); b'^l^ = ^{^B,{kpq) + B^ikpq)); 

i 

bZ, = i(^i(^w) + B2{kpq)), (A12) 

and substituting into Eqn. (A9), we obtain immediately the result of Eqn. (11). 

Finally, we derive the energy equation in the moving reference frame. Using conservation of energy 
triad by triad (which is preserved, together with incompressibility, when passing to the Lagrangean reference 
frame), we can write: 

DtCu{t)\t=o ^^ I dpdqkpqekp,[-Hikpq) + H{pkq) + H{qkp)l (A13) 

where: 

H{kpq) = [6fcp, + Wj,b'Z^]CkC, ^ w.b^^CpC,. (A14) 

We obtain immediately: 

DtCk{t)\t=a ^ J~ j (ip<iqkpqekpq[{bkpq + WqbZ^)CpCq 

-ibkpq + {^rWqb'^;^^^^)CkCq], (A14) 

where: 

K'L = ^(^B2{kpq) + B,{kpq) + Bi{kpq)). (A15) 

A form of Eqn. (A14), which is manifestly zero at cquipartition, can be obtained using the relations: 

B2{kpq) = B2{p,k,-q) and Bi{kpq) = Bi{p,k,-q). (A16) 

Substituting Eqns. (A15-16) into Eqn. (A14), we obtain finally the result shown in Eqn. (15). 

Appendix B. Noise correlation 

In order to calculate the amplitude of the energy flux fluctuation gn {m) , we must evaluate the contrac- 
tions of the product: 

i;"(zi)K(zi + n) - «''(zi)]9;3«"(zi + ri)t;^(z2)K(z2 + ra) - v%^2)]dav'^ {^2 + rs), (Bl) 
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which are Hsted in Table Bl. In this section, we shall adopt the notation: V'a = dai^- Of these contractions, 
numbers (Bl.l), (B1.4), (B1.7), (Bl.lO) and (B1.13) do not contribute because of incompressibility. The 
remaining ones fall into two groups: (B1.2), (B1.3), (B1.4) and (B1.6) are in the form: < g{0)g{t) >^< 
t'(0)w(0) >< v(0)v{t) >< v{t)v{t) >, where < v{0)v{t) > is associated with velocity fluctuations on the scale 
of the wavepacket radius R. These terms come from noise components of g in the form v < (vV)v >, which are 
identicaly zero, and (wV) < w^/2 >, which are associated with sweep and do not contribute to energy transfer. 
The second term contains all the other terms, which are in the form < v{0)v{t) >< v{0)v(t) >< v{0)v(t) > 
and act on the same scale of fci and fc2 . 

It is possible to show explicitly that the contractions in the first group do not contribute to < g^ >. 
Consider for example contraction (B1.3): 

< g„i (mi, ii)5„2 ('7^2,^2) >''^'= / d^rid^r2H{ri,ni)H{r2,n2) 
X / d^zid^Z2Wizi,mi)W{z2,m2)[C°'^iri,0) - wC^iO, 0)] 



: [C7(zi + ri - Z2 - r2,t) - wC^^z^ + n - Z2, i)]CJ^(r2, 0) 

(tr (S^ ^l]^ C,C,C,eM~Vr>\t\)W,{rn^)W,{rn2) 
f (kp)2(pq) (kp)(kq) -( 

I /c2p2 fc2 / 

x[(2 - Wk{k+p))H{k+p- fc„J - Wk{k+p)H{p - fc„J] 

X [(2 - Wp{p + q))H{p + q- fc„J - Wpik)H{q - fc„J]. (B2) 

The geometric term in braces is antisimmetric in q, so that that only nonzero contributions are proportional 
to H{p + q — fcna)- We change first variables: q + p ^ q, so that H{p + q — kn^) — > H{q — fcnj). Because of 
the factor Wp{mi)Wp{m2) , we can expand in a power series in p: 



-k!^dk^ 



(M-^/^f^^^C%^-»^)l -4(M-V(...-.V) (B3) 



V k^p'^ k^ / 



3 



Indicating by u the cosine of the angle between the planes kp and kq, we have: 



xyz — X z — uxzy'{l — x2)(l — z^), 

so that, substituting Eqn. (B3) into (B2) gives zero. Going to next order produces a factor p^x^, which 
gives zero again, while still the next order leads to a contribution which is by a factor (fci?)^" smaller than 
the terms in the second group. 

The calculation of the terms in the secon group is very tedious but straightforward and follows the same 
lines of Eqns. (28) and (B2). We conclude by giving the exact form of the coefficients entering the expression 
for the fluctuation amplitude < g^ > [Eqn. (30)]: 

Bi{kpq) = —pkyz {y + xz); B2{kpq) = kp {1 + z ){z + xy); 

B^{kpq) = kpxy (1 - z^). (B4) 

Hi = [(2 - Wp{k))H{k - fc„J - H{q - k^,)wp{k)] 

X [[(2 - Wq{k))H{k - fc„J - H{p - kn,)Wq{k)] 

-[(2 - w,{p))H{p - fc„J - H{k - k^,)w,{p)]]- (B5) 

H2 - [(2 - Wp{k))H{k - fc„J - H{q - knJWpik)] 

X [[(2 - Wp{k))H{k - fc„J - H{q - kn,)wp{k)] 

-[(2 - Wp{q))H{q -ka,)- H{k - k^,)wp{q)]]; (B6) 

H^ - [(2 - Wp{k))H{k - fc„J - H{q - k^,)wp{k)] 
x[[{2-Wk{q))H{q- fc„J - H{p - kn^)wk{q)] 
-[(2 - Wk{p))H{p -ka,)- H{q - ka,)wk{p)]]. (B7) 

21 



Appendix C. Coefficients for the shell energy equation 

The RHS of Eqn. (52) contains the following contributions: 

A term giving the effect of the timescalc becoming shorter at large n: 



-F At I dn g{n,n' ,t) F{n')dl,E{n' ,t) 



where: 



F dt I dn'g{n,n',t){l-e) < e„(0)e„+„'(t) > 



S(n',t) = ^i4^^|^±4^, and F(n') =< CnWCn+nKO) > 



<e„(oK„+„'(o)> 

A term coming purely from the noise, due to Beta Model effect: 



(CI) 
(C2) 



F dt dn'g(n,n'-n,t)<(C„(n,0)-C„(n-l,0))(e„.(n',i)-C„.(n'-l,t))> 



F(3 dt I dn'g{n,n',t) < ^„(0)C„+„-(t) >; 



where: 



13 = 1 + a-^ -2^(1 + 0^)-^. 
A cross correlation term between Beta Model effect and timescale mismatch: 



(C3) 

(C4) 



2F dt dn'g{n, n' , t)d'„E{n' , t) < [e„(n, 0) - in{n - 1, 0)]C„+„'-i(0) > 



F(i' / dt / dn'g{n,n',t) < e«(0)e„+n'-i(i) >, 



where: 



/3' = 27 1- 



1 + a 



-2^, 3 



(C5) 
(C6) 



Finally there is the correction to the relaxation terms of the shell equation, from Beta Model effect: 
(■D--)< ((/)„(n,O)-</)„(n-l,O))(</)„(n,i)-0„(n-l,t))> 

To lowest order, there are no cross correlation between {D — ■j)(0g'("-) ~ 't'oin + 1)) and previous terms due 
to the origin of (J)g from 9n'C rather than from ^. 

Using the formula J^ dt exp(— ^ — ct^) = y^exp(— 2v6c) [36], and shifting the n' integration so that 
n ^ 0, the integral term in Eqn. (52) can be rewritten in the form: 



(C7) 



dn' 



d^ 



<eo(0)e„'(0)> [7'(^-l)+/? 



+f3' <^omn'-l{0)>jr- 



4c2 



exp((i — 2v5c), 



(C7) 



where: 



/ttD 



b-—- 
~ AD' 
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c= — +cxp(7n); 



d^ 



n'V 
'2D' 



(C8) 



Writing explicitly: 

|dn'{[7^( 



1 - 



32D'^ 



6D\n'\ 



n'2 



■/3 



(T/2 + 415)2 (v-2 + 4£,) i y2 + 4£, y ' ^ J 1 + 50(a _ i)2„/2 
,,4D+\n'W'iD + V^ 



+P' 



AD + V^ 



I I 

l + 50(a-l)2(n'-l)2J 



F 



V4JDTV^ 



exp 



2D 



(n'- \n'\y/4:D + V^) 



(C9) 



It remains to evaluate the term D < (j>A'^4' > on the LHS of Eqn. (51): 

< 0A2</> >== / dtidia / dnidna F{ni - n2)g{ni,ti)dl^g{n2, ^2) ^[h - ^2] 

= 51 ~T / '^W dnidn2F(ni - n2)g{ni,ti)d^dl^g{n2,t2), 



where 



/"dii™S(i). 



(CIO) 



(Cll) 



Using the defining relation for g: {dt — Dd^-^ + Vdn)g{n,t) ~ S{n)S{t) (in the continmim limit) and the 
condition: g{n, 0) = g{n, cx)) — 0, we can write: 



00 „ ^00 
ml ./o 



m—0 



< (/.A20 >^Y,^ I dt / dnidn2 i^(ni - n2)g{nuh){Ddl^ - Vdn,rdl^g{n2, ^2) 



00 ,— , 






2L>(nf + n|) 



2N2- 



Integrating by part we obtain a series in the form: 



< 0A20 >= J2 



m=0 



dnidn2 



2D{n\+nl) 



{a^mdl^ + hradl'^-^)F{ni - n2). 



1\2 V"'™"n2 



(C12) 



(C13) 



Given the form of the correlation F [see Eqn. (43)], the action on it of the operators (amd^™ + bmdf™ ^) 
will lead to terms in the form: 

Ci+C2(Tli-n2) 

(l + 50(a-l)2(ni-n2)2)p' ^ ' 

Substituting back into Eqn. (C13), we see that terms in ci do not contribute because of the antisimmetry 
of the integrand under the transformation ni —> n2, n2 -^ rii. Similarly for the terms in C2, due to 
the antisimmetry under simultaneous change of sign of ni and 712. Hence the term < (/)A2(/) > does not 
contribute to Eqn. (51). 
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FIGURE CAPTIONS 

Fig. 1. Kolmogorov constant CkoI and dimensionless parameter p, as measured from the energy flux through 
scale k, in a reference frame moving with velocity v(A', k) — w(A', k)v, < A' < 0.9, with A = 0.9, fixed 
[ see Eqn. (14)]. 

Fig. 2. Energy transfer profiles T(k,p) for fc = 1. (a): Lagrangean frame using Eqn. (15); A' — 0.9. (5): 
laboratory reference frame, using Eqn. (10). (c): Lagrangean frame using simplified closure with 
Eqn. (17) and: Wp^q = Wp,q{0.7,k). (d): the same as (c), but with Wp — Wp(0.7,imii(k,q)) and 
Wq = Wp{0.7,m.m{k,p)). 

Fig. 3. Normalized noise correlation — " "at^ — vs. shell separation x = {a — l)n'. (a): calculation from Eqn. 
(30); (b): fit of (a) by: F{x) = ("l + SOx^)-!. (c): profile for A{x) = J dp\T{k,p)T{k + x,p + x)\, 
normalized to A{0); 

Fig. 4. Finite difference coefficients Ar vs. a [ see Eqns. (33), (36) and (37)]. Notice that higher order differences 
are always dominated by the advection-diffusion part. 

Fig. 5. Comparison of gaussian part of the fluctuation amplitude << vf >|> — << vf >^>, obtained from 
integration of Eqn. (49) (a), and from direct average of w*"' (b). 

Fig. 6. Scaling correction to S'(Z,a, c^,4), as a function of a for c^ ~ 2; notice the saturation occurring at 
a~ 1.3. 

Table Bl. Contractions contributing to < gnijn)gn'(rn!^ >. 
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